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Abstract 

We consider risk-neutral returns and show liow tlieir tail asymptotics 
translate directly to asymptotics of the implied volatility smile, thereby 
sharpening Roger Lee's celebrated moment formula. The theory of regular 
variation provides the ideal mathematical framework to formulate and 
prove such results. The practical value of our formulae comes from the 
vast literature on tail asymptotics and our conditions are often seen to be 
true by simple inspection of known results. 

1 Introduction 

Consider risk-neutral returns X with cumniulative distribution function F . We 
will impose a mild integrability condition on the right and left tail denoted by 
(IR), (IL) respectively. We write F = 1 — F and, if it exists, / for the proba- 
bility density function of X . The class of regularly varying functions at +oo of 
index a is denoted by R^. The reader not familiar with the theory of regular 
variation should think of positive multiples of and harmless perturbations 
such as (l -I- x"/^) , a;" log a; or x"/ (logx)'^. By convention, functions in Ra 
are (eventually) positive. 

The purpose of this paper is to connect the tail behaviour of X to the wing 
behaviour of the Black-Scholes implied volatility, sharpening Roger Lee's cele- 
brated moment formula [9] . From a mathematical point of view, the challenge 
is to relate the asymptotics of the distribution F to the asymptotics of a non- 
linear transform, namely the Black-Scholes implied volatility. From a financial 
point of view, we give further justification of volatility smile parametrizations 
seen in the industry and obtain new insights into the wing behaviour of a variety 
of models. 



The normalized price of a Black-Scholes call with log-strike k is given by 

CBS (fc,C7) = $(di)-e^$(d2) 

with di_2 {k) = —k/a± a/2. The implied volatility is the (unique) value V {k) 
so that 

/•OO 

CBS {k, V (fc)) = / (e^ - e") dF (x) =: c (fc) . 

J k 

We define the strictly decreasing function ij) : [0, oo] — > [0, 2] by 

■)p[x] — 2 — 4 \/ + X — X . 
Theorem 1 (Right-tail- wing formula) Assume a > and 

3e > : E[e(^+')^] < oo. (IR) 

Then 

(i) => (ii) => (in) => (iv), 

where 

-log/(fc)Gi?„; (i) 
-log^(fc) e (ii) 
-logc(fc) e -Rai (iii) 

and^ 

V{kf/k ~ V [- log c (fc) /k] . (iv) 

// (ii) holds then — \ogc{k) ^ —k — logF and 

V{kf /kr^ip [-1- log F{k)/k], (iv') 
if (i) holds, then — log / — logF and 

F(A;)Vfc~ V'["l-log/(fc)/fc]. (iv") 

Finally, if either — log / (fc) jk or — \ogF (fc) jk or — logc (fc) /fc (?oes to infinity 
as k ^ oo then V"^ {k) behaves sublinearly. More precisely, 

2 1 1 1 

^^^^ " -21og/(fc) /fc"'' -21ogF(fc) /fc -21ogc(A;) /fc' 

We emphasize that (iv),(iv'),(iv") contain the full asymptotics of the implied 
volatility smile. For instance, we can see when lim sup V (fc)^ /k in Lee's moment 
formula is a genuine limit: V {k) /k converges if and only if —\ogF{k) jk 
converges to some limit Q. Note that our condition (IR) forces ^ > 1. Note also 
that in this case — log F G Ri so that condition (ii) is automatically satisfied. 

(fc) ~ h (fc) means g (fe) /h (fe) — » las k — » oo. 
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In models without moment explosion (Black-Scholes, Merton's jump diffu- 
sion model, FMLS with (3 = —1...) the moment formula indicates sublinear 
behaviour of the implied variance, limsup V{kY /k = \miV{kY /k = 0, but 
yields no further information. In contrast, theorem 1 gives the precise sublinear 
asymptotics. For instance, in the Black-Scholes model log Jbs {k) ~ —k'^/ (2o-^) 
and (v) implies V (k)"^ ~ cr^ in trivial agreement with the Black-Scholes flat 
volatility smile. 

There is a similar result which, as fc ^ oo, links / (— fc) , F {—k), normalized 
out-of-the-money put prices p{—k) and the implied volatility in the left wing. 

Theorem 2 (Left-tail-wing formula) Assume a > and 

3e > : E[e-'^] < oo. (IL) 



Then 
where 



and 



(i) => (a) => (Hi) => (iv), 

-log/(-fc) e (i) 

-logF(^fc)ei?,; (ii) 

-\ogp{-k) e Ra, (iii) 

F(-A;)Vfc~ V[-l-logp(-fc)/A;]. (iv) 



// (ii) holds then — logp (— fc) ~ — logF (— fc) and 

V{-kf/kr^ip[-logF{-k)/k], (iv') 
if (i) holds, then — log / (— fc) — logF (— fc) and 

V{kf/k^i,[-\ogf{-k)/k]. (iv") 

Finally, if either — log/(— fc) /k or — logF{—k) /k or — \ogp{—k) jk goes to 
infinity as k ^ oo then {—k) behaves sublinearly. More precisely, 

V{-kf/k ~ _2iogf{-k)/k°^ -21ogF(-fc) /k '^^ -2\ogp{-k)/k' ^^"^ 

In conclusion, under mild integrability and regular variation conditions, tail 

asymptotics translate directly to asymptotics of the implied volatility smile. The 
practical value of formulae (iv'), (iv") comes from the vast literature on tail 
asymptotics and our conditions are often seen to be true by simple inspection of 
known results. The authors would like to thank Jim Gatheral, Roger Lee and 
Chris Rogers for related discussions. Financial support from the Cambridge 
Endowment for Research in Finance (CERF) is gratefully acknowledged. 
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2 Elements of Regular Variation Theory 

Definition 3 A positive real-valued measurable function f is regularly varying 
with index a, in symbols g & Ra if 

lim ^ = 

x^oo g[x) 

Functions in Ro are called slowly varying. 

The definition extends immediately to functions whieh are positive for x large 
enough. Recall that g ^ h means lirsix^oo g{x)/h (x) = 1. Trivially, g G Ra and 
g implies h G Ra- The following result can be found in the fine monograph 
[5, Thm 4.12.10, p255]. 

Theorem 4 (Bingham's Lemma) Let g € Ra with a > such that that 
is locally integrable at +oo. Then 

/■oo 

-log / e-^^^Uy g {x) . 

J X 

3 Right- Wing Smile Asympotics 

Proof of Theorem 1. We first remark in presence of condition (IR), under 
either assumption (ii) or (iii) one must actually have a > 1. Indeed, 

F (k) = F[X >k]< e-(i+^''=E[e(i+^'^] =^ - logF (fc) > (1 + e) k. 

Assuming — logF e i?a for a < 1 would entail — log F{k) < k"' for a' £ 
(a, 1) and k large enough which contradicts the lower bound just established. 
Similarly, 

c(fc) = E[(e^ - e'=)+] < E[e^;X> k] < E[e^e<^-''^] = e-'''E[e^^+'^^] 

so that — logc (A;) > ek and — logc G Ra can only happen with a > 1. 

(i) =^ (ii) : Use Binghan's lemma with the (eventually) positive function g (k) = 

— log/ (A:) e Ra- (ii) (iii): Lemma 6. (iii) => (iv): Lemma 7. The final 

statement on sublinear asymptotics follows from ip[x] 1/ (2a;) as a; — > oo. 

Indeed, 

^[.]=2-4.(yrTT7^-l)=2-4.(^-^+0(.-3))^l. 



Lemma 5 Assume (IR)- Then the normalized call price with log-strike k is 
given by 

f-OO 

c{k)= / e^F(a;)dx. (1) 

Jk 
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Proof. Note that F = 0(6"^) since F (x) = F[X > x] < e-(i+^)^E[e(i+'^)^]. 
Integration by parts gives 



c{k) 



POO 

7 

Jk 



') dF (a 



= -f'(x)(e^-e 
and the boundary contribution disappears. 



k\ 1 00 

k 



e'^F (x) dx 



Lemma 6 Assume (IR) and that — logF e Ra for some a> 1. Then — logc G 
Ra and as fc —> 00, 

-logc(fc) ~ -fc-logF(A;). (2) 

Proof. Obviously k ^ k £ Ri and — logF G Ra,c( > 1. We want to apply 
theorem 4 with k^ ip{k) = — logF {k) — k to obtain (2) but need to be careful 
since in general the difference of two regularly varying functions may not be 
regularly varying. From the last proof we know that — logF(fc) > (1 + e) fc. 
When Q > 1 then (f (k) is immediately seen to be in R^- When a = 1 we can 
write — log F {k) = kL (k) with slowly varying L{k) > 1 + e. To see that tp G R^ 
it suffices to write 

tp{Xk) L{\k)-l L{\k)/L{k)-l 



1 



so that 



ip{\k) 


< 







^{k) L{k)-1 

L (Afc) /L {k) - 1 



1/L{k) 



< ^-^\L{Xk)/L{k)-l\ 



l-l/L (fc) 

and this tends to zero as fc ^ +00 since L € Rq. In either case, we can apply 
theorem 4 with g = if and obtain 



- log c (fc) = - log / e^+'°s ^(^)da; ~ - (fc + log F (A;)) = 

Jk 



(fc). 



Note that (p (•) was seen to be in R^ so that — log c (fc) ~ (k) must also be in 
Ra- ■ 

Lemma 7 Assume (IR) and — logc e Ra for some a>l. Then, as k ^ 00, 



log c (k) 
k 



and 



k 1 _ V{k)^ 
'' ~'2V{k)'^ ^ 2^ Sk 

V{kf /-logc(fc) 



O 



logfc 



k 



k 



Proof. Appendix. 



Remark 8 In the preceding lemma the condition of regular variation can he 
replaced by the weaker 



3n e N : lim inf V{k)k"- > 1. 

/>;—>■ 00 
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4 Left- Wing Smile Asympotics 



The proofs are similar and therefore omitted. 

Lemma 9 Assume (IL). Then the normalized put price with log-strike —k is 
given by 

/—k POO 
e''F{x)dx= / e-''F{-x)dx. 
-co J k 

Lemma 10 Assume (IL) and that — logF(— fc) is regularly varying. Then 
— logp {—k) is regularly varying as k ^ oo and 

-logp{-k) k-logF{-k). (3) 

Lemma 11 Assume (IL) and that — logp(— fc) is regularly varying. Then, as 
A; — > oo, 

k \ k 

5 Examples 

In practice, one has X — log {St /Ft) where St denotes the risk-neutral stock 
price at time T and Ft is the time-T forward price. Also, k = log (K/Ft) and 
all quanitities f,F,c,p,V depend on time T and we set V (k) = V {k,T) =: 
f^BS (k, T) -\/T, now calling obs {k, T) the implied volatility. It is worthwhile to 
spell of part (iv) of Theorem 1, 

a%s (fc,T)T/fc~V[-logc(fc,T)/fc]. 

Noting that {St /Ft ■ T > 0) is a martingale and using convexity of the call 
payoff, it is easy to see, for k fixed, c (fc, T) —E (e^ — e*^) is non-decreasing in 
T and so is ip [— logc(fc,T) /fc]. In fact, J. Gatheral and E. Reiner^ point out 
independently that the implied total variance a%g (fc, T) T is non-decreasing in 
T, a consequence of monotonicity of the undiscounted Black-Scholes prices in 
a^T. Thus, our asymptotic results respect the term structure of the volatility 
surface. 

In the examples below, we will focus mainly on applications of the right-tail-wing 
formula, applications of the left-tail-wing formula being nearly identical. 

5.1 Sanity Check: Black-Scholes Model 

If (T denotes the Black Scholes volatility, the returns have a normal density with 
variance cr^T. Obviously then, 

log /bs (fc)~-fcV (2<7'I^) 

and Theorem 1 implies ass {k, T)"^ ^ as fc — > oo in trivial agreement V = a. 
^Presentations at the Global Derivatives & Risk Management Conference 2004. 
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5.2 BarndorfF-Nielsen's NIG Model 



Here X = Xt ^ NIG (a, /3, /uT, 5T). The moment generating function is given 
by 

M {z) = exp T (^5 I ^o?-0^ - ^a^ - (/3 + 2)^| + /z^ 

It is custom to write a = -\//?^ + 7^ with 7 > 0. Prom [4] and the references 
therein we have 



/ (fc) ~ C |A;|-'/' e-V^^|fc|+/3fe as ^ ±00 
and we see that — log / is regularly varying (with index 1). Moreover, 



log/(fc)/fc^ 
and from Theorem 1 

k 



/3 + 7^ + /3 as A: ^ +00. 



^^(-l-log/(fc)/fc) 



r+72-/? 



in agreement with Lee's moment formula with critical moment 1+p* = \/ (3^ +7^ — 
/3 — 1 = a — /?, as can be seen directly from the moment generating function, 
see [9]. 

5.3 Carr-Wu's Finite Moment Logstable Model 

Here X = Xt ~ La (^r, oT^I'^ ^ — l) where the law La {0, cr, (3) has characteris- 
tic function 

|g j^^mX j ^iu9—\u\'^ (7^ (l — «/3(sig u) tan ■^). 

with a e (l,2],6l e M,fT > 0,/3 e [-1,1]. From [6, page 10, equation (6)] and 
the references therein, 

- logi?' (fc) ~ X [Taa" |sec (7ra/2)|]"^/^""^' as fc ^ 00, (4) 

and from Theorem 1 we see that 



als{k,T)T ^k^-^ X ^[Taa"|sec(7ra/2)|]^/("-^^ ^Oasfc 



00. 



Note that in the limit a1 2 the Black Scholes result is recovered. We also note 
that the moment generating function of La {6,a,—l) exists for all positive z 
and is given by 

M(z) = exp z^^ — (zct)" sec 
Kasahara's exponential Tauberian theorem [5, p253] then gives immediately (4). 
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5.4 Merton's Jump Diffusion Model 

As in the examples above, the return process X. is Levy with triplet (/i, cr^, K) 
where if is A (=intensity of jump) times a Gaussian measure with mean a and 
standard deviation S describing the distribution of jumps. We now demonstrate 
how to proceed without explicit knowledge of the asymptotic tail. Set X = Xt 
and note that E [exp [zX)] < oo for all z. Optimizing over z we get the tail 
estimate 

F{k) = F[X > k]< inf e-^*=E [exp (zX)] = e^(^'")-^*'= 

z 

where K (z) = log E [exp (zX)] is the logarithmic mgf and z* = z* {k) is deter- 
mined from^ 

K' {z*)=k. 

For the Merton model, 

K{z) = T l^/x + + A (e^«+-'^V2 _ I 

from which for ^ > it is easy to see that 

\/21ogfc 

z =z {k) 

so that K (z*) - z*k ^ -z*k and 

k 



logF(fc) < -z*k = --^/2\ogk. 

From the saddle point results of [7] or the Levy tail estimates from [2] specialized 
to this example, 

logi^(fc)~-^V21ogfc, 

and Theorem 1 implies 

ala(k,T)T^6 X , ^ 

Note that this is independent of the mean jump size a provided S > 0. If 5 = 
and a > a similar argument shows that z* = z* {k) logfc/a and 

4,(fc,T)T~ax^. (5) 

Remark 12 (J. Gatheral) A Poisson process with intensity A has n jumps 
with probability e~^\^/n\ = e^^"^ with g(n) ~ — nlogn by Stirling's form,ula. 
The time T = 1 Black-Scholes value of a digital is $(^2) and log$(cJ2) ~ 
— fc^/2cr^. Identifying k na leads to 

k 

-k^ fta^ ~ - (klo) log A; =4- ~ a x — 

21ogfc 

in agreement with (5). 

^Readers familiar with large deviation theory recognize the Fenchel-Legendre transform. 
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6 Further Examples and Discussion 



Our tail wing formulae apply as soon as one has some (in fact: crude) asymptotic 
knowledge of the tail behaviour of the retTirns. As remarked in the introduction, 
for many models fine tail asymptotics are available in the literature and the smile 
asymptotics follow. In particular, the estimates of Albin-Bengtsson [1, 2] allow 
to use our results for the vast majority of exponential Levy models [10]: NIG 
appears as a special case of generalized hyperbolic processes, Meixner processes 
as special cases of GZ processes. CGMY and Variance Gamma model are also 
covered. Tail estimates for stochastic volatility models appear in the literature, 
although one has to be careful with results obtained in short time regimes since 
the limits T ^ and fc — ^ oo in general do not interchange^. We remark that 
condition (IR) rules out models for which every p*''-moment of the underlying, 
p > I, explodes and a similar remark applies to condition (IL). Unfortunately, 
there are stochastic volatility models with this degenerate behaviour [3] and our 
(current) results do not apply. 

When a moment generating function is known, one can often use an exponential 
Tauberian theorem to obtain log-tail estimates, as demonstrated in the Example 
5.3. Also, the Fenchel-Legendre transform gives quick upper bounds which often 
turn out to be sharp, compare Example 5.4. We finally remark that some of 
the implications in Theorem 1 and Theorem 2 can be reversed under Tauberian 
conditions but we shall not pursue this here. 



7 Appendix 

Proof of Lemma 7. The following bounds for the normal distribution func- 
tion $ are well-known and can be obtained by integration by parts (or other 
methods), 

-(1 - ^) < H-x) < ^7^, x>0. 



/2ttx V27ra; 
Using assumption (IR), 

E[{e^ - e'=)+] < E[e^;X > k] < E[e^e<^-''^] = e-''''E[e^^+'^^] 

we see that c (fc) is exponentially small in k, and this implies that 



Vlk)"^ 

a:=limsup ^ ' <2 (6) 



In particular, both 



k V{k) k V{k) 

di(k) = — H —, d2(k) = — — 

^ ' V{k) 2 ' ^ ' V(k) 2 



^The wrong wing behaviour in Hagan's SABR formula, based on short time asymptotics, 
is a warning example. 
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will then be negative for k large enough so that we will be able to use the bounds 
on From the Black-Scholes formula for a normalized call with log-strike k, 



c(fc) = $(di)-e'=$(d2), 



we obtain 



and 



27rc{k)<e-</^-j- 



'^e-'^l/^ U±] 1-']. 



Because df/2 = d^/2 — k, this simplifies to 



-4/2 



1 



dl 



+ 



do 



We now define ei = ei (fc) by 



< V27rc {k) < e 



dl 



-dl/2 



4 



_ 1 1 

dl d2 



4 



logc(fc) = -^+ei(fc) 



and will show |ei {k)\ = O (logfc). To start, we note the bounds 

dl df} d2 ~ - di'^ d2 dl 

From (6) we know that V{k) < \/2fc, hence 



-d2 = 



k 



V{k) 



V{k) 2 
Therefore, for k large enough, 



> \fkl2 — > oo as fc — > oo. 



(7) 



^e^m < _ J_. 

dl 

In fact, (6) gives a! G (0, 2) s.t. V (fc)^ < a!k for all k large enough. From the 
definition of di we then see that 



a"'\/fc := ^ 



1 



< -dl {k) 



with a'" > because a' G (0, 2). Hence 

2^e'^('=) < 
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and wc sec that ei (fc) — > — oo as fc — > oo and we only need a bound on — ei (fc). 
We start by showing 

3n e N : lim inf V{k)k'' > 1. (8) 

/C— '■OO 

To see this note that for k large enough ei < so that 

-logc(A:) = |-ei(fc)>| 

By assumption — log c G i?a for a > 1 so there exists L G Rqso that — log c (fc) = 
k^L (k) and 



so that 

< 2k°'L {k) + k = k'^L (fc) 



where L G Rq and ^ (fc)^ /c" ^ = 1/L (/c) and (8) holds with any integer n > 
(a -2)+. 

Wc return to establish a bound on — ei (fc). We already have lower and upper 
bounds on the implied volatility, namely (replace n by n + 1 if needed) 

^-2n ^ Y j-^-j2 ^ ^.^j^ a' <2 and /c large enough, 
which we can use to bound di^2-, namely 

-dl < A;^+" - ^ < (9) 



-d2 < fci+" + ^<2fci+" 
2ft; 

(at least for A; > 1). In order to derive an upper bound for — ei recall that a 
lower bound on \/27re^^ was given in (7) by 



£2 (A;) 



rflo dl 





















1 

+ -1^ 



(_fc + y2/2)%(_fc_y2/2) _ 

We already noted the existence of a' € (0,2) such that := < (3k < k 

with (5 = 12 < 1 for fc large enough. Noting that, for k fixed, the function 

[0, 9 a; {-k + xf + (-A; - x) 



11 



is strictly decreasing on [0, (3k] we can get a lower bound on 62 as follows, 

{-k + f3kf + {-k-pk) 



(*) > 



dld2V 
dld2V 

> ^=—5 (lor k large enough) 



> 



2 



= ^fc2-4(i+")-i/2 with 7 > 0. 

Therefore, 

\/2^e^i('=) > 62 (fc) > 7fc-^/2-*" > 
and after taking logarithms 

-ei (fc) < (4n + 5/2) log A; + C 

for some constant C = C (7). In particular, 

|ei(fc)| =0(logfc). 

Recall that ei was defined by 

logc (fc) = -y + ei(fc) = ^ - ^ - + ''^^^ 
and after dividing by k we have, 

logc(fc) /I fc 



(^(logfc/fc) 



2 



k \2 21/2 8/0 
which tends to zero as fc — > 00 so that 

logc(fc) 1_ k _V 
k 2 ~ 2^2 ~ 8fc' 

We can rearrange the above equations as a quadratic equation in V^/k, 

logc(fc) ei(fc) _ 1 V^/k 1 
k ^ k ~ 2y2/fc + 8 2' 

One solution is given by 

Yl=^ /_logc(fc) ^ ei(fc)\ ^ ^ / logc(fc) 
/l y k J \ k 

while the other is rejected as it violates (6). ■ 
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